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Qh 1 We calculate the complete next-to-leading order (NLO) QCD corrections to the J/ip + W as- 

sociated production within the factorization formalism of nonrelativistic QCD at the 7 TeV LHC. 
We provide the numerical results for the leading-order (LO), NLO QCD corrected differential cross 
sections of the J/ip transverse momentum by adopting the event selection criteria requested by the 
ATLAS experiment. We find that the differential cross section at the LO is significantly enhanced 
by the NLO QCD corrections. 
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Recently, after we published our work on the QCD corrections to J/ip production in association 
with a H^-boson at the 14 TeV LHC [I], Darren Price in the ATLAS experiment suggested us to 
extend our study to the colliding energy of 7 TeV, and asked for theoretical predictions for the 
corresponding differential cross sections in the restricted kinematic regions of the final particles. In 
Ref.jTj we presented the detailed calculations of the NLO QCD corrections to the J/ip + W associated 
production at the -y/s = 14 TeV LHC, but not for the 7 TeV LHC. In this letter, we calculate the 
NLO QCD corrections to the associated J/ip production with a W gauge boson in the nonrelativistic 
QCD (NRQCD) at the 7 TeV LHC, and provide the theoretical predictions for the px distribution of 
J/ip by adopting the event selection criteria requested by the ATLAS experiment. 

We calculate the pp — > J/ip + W + X process by applying the covariant-projector method [2] in the 
NRQCD framework. Some representative Feynman diagrams for the partonic process ud — > J/ip + W + 
at the LO and the QCD NLO are depicted in FigfTJ The cross section for the pp —> J/ip + W + X 
process is expressed as 

Ef (0 H \n\) 

/ dx 1 dx 2 . T \ ^_ J/ , . o- (i + j -»■ cc[n] + W, s = X!X 2 s) 
,.., J N col (n)N pol (n) 

x [G i/A (x 1 ,n f )G j/B {x 2 ,fif) + {A^ B)] , (0.1) 

where (ij) = (ud) and (ud), N co [(n) and N po i(n) refer to the numbers of colors and polarization of 
cc[n] states produced [2j, respectively, a is the cross section for the short distance production of a 
cc quark pair in the color, spin and angular momentum state n at the parton level. (0 J '^[n\) is 
the long distance matrix element, which describes the hadronization of the cc quark pair into the 
observable charmonium state J/ip. Gij/ab are the parton distribution functions (PDFs), A and 
B refer to incoming protons at the LHC. We ignore the processes which are suppressed by a tiny 
Kobayashi-Maskawa factor. 

In the leading order (LO) calculations, only the Fock state 3 S\ is involved. However, for the 
NLO QCD corrections to this process, we need to consider both the virtual corrections and the real 
gluon/light-quark emission processes. The virtual corrections only contain the contribution from 
cc Fock state 3 S\ , while the real gluon/light-quark emission processes should involve 1 Sq , 3 5| 
and 3 Pj Fock states. In our calculations, we apply the dimensional regularization (DR) scheme to 
regularize the UV and IR divergences, and the modified minimal subtraction (MS) and on-mass-shell 




Figure 1: Representative Feynman diagrams for the partonic process ud 
the QCD NLO. 



J/xp + W + at the LO and 



schemes to renormalize the strong coupling constant and the quark wave functions, respectively. In 
FigEJ we present the divergence structure and divergence cancellation routes in the NLO calculation 
for the pp — > J/ip + W + + X process. In the virtual correction calculations, we find that only Figfljf ) 
and FigQTh) induce Coulomb singularities, and we use a small relative velocity v between c and c to 
regularize them [3]. The two cutoff phase space slicing method (TCPSS)[4J has been employed for 
dealing with the soft and collinear singularities in real gluon/light-quark emission corrections. 

We take CTEQ6L1 PDFs with the one-loop running a s in the LO calculations and CTEQ6M PDFs 
with the two- loop a s in the NLO calculations [5], respectively. For simplicity we define fj, = fi r = [if 
and take the renormalization/factorization and NRQCD scales as \x = m T or myj, fi\ = m c , 



respectively, where m T = \l [p T ) + rn 'j/^ is the J/ip transverse mass. The masses of the external 



particles and the fine structure constant are taken as raw = 80.398 GeV, m. 



c — 2 nl J/ip 



1.5 GeV 



and a = 1/137.036. We use two sets of color-octet (CO) long distance matrix elements (LDME) listed 
in Table [T] as the input parameters, which were obtained by fitting the differential cross section and 
polarization of inclusive J/ip production data from various experiments [6j[7]. 

We calculate the LO and NLO QCD corrected cross sections for the pp — > J/t/} + W — > J/ip + 
M + v p, + ^ process at the 7 TeV LHC by adopting the following kinematic cuts on the final particles 
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Figure 2: The divergence structure and divergence cancellation routes in the NLO calculations for the 

pp — > J/ip + W + + X process. 





Set 1 (GeV' A ) 


Set 2 {GeV' 6 ) 


< o J ^[ 3 s[ s) ] > 

< 0-WppW] > 


3.04 x 10~ 2 

1.68 x 10~ 3 

-9.08 x 10" 3 


8.9 x 10~ 2 

3 x 10~ 3 

5.6 x 10~ 3 



Table 1: Two sets of CO LDME for J/ip. Set 1 are extracted via a global fit from various hadroproduc- 
tion, photoproduction, two-photon scattering and electron-positron annihilation experiments [6j. Set 
2 are extracted by fitting the differential cross section and polarization of prompt J/ip simultaneously 
at the Tevatron [7]. 



according to the requests of the ATLAS data analyses: 



p^ > 25 GeV, \if\ < 2.4, p u T > 20 GeV, mf > 40 GeV. 



(0.2) 



where mW is W boson transverse mass defined as 
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^PtPt [1 — cos l 



1, and (f)^ 



is 



the angle between muon and neutrino in the transverse plane [8]. 

In our calculations the uncertainty for theoretical prediction mainly comes from two parts: the 
short distance part and the long distance part. For the former part, the dependence of the cross 
section on the renormalization scale [x T and factorization scale fif induces theoretical uncertainty. 
For the later part, the values of CO LDME for the J/i/j extracted from the experiments by different 
experimental groups vary significantly. In the following tables and figures, we will give the results by 
adopting two set of CO LDME and two typical scales (/u = m T , mw)- 

In Table [2] and Table El we provide the cross sections for the pp — > J/ip + W — »■ J/ip + fi + v^ + X 
process at the 7 TeV LHC in the p T and the J/ip rapidity intervals by taking the two different 
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Table 2: The cross sections for the pp — > J/ip + W — > J/ip + \x + v M + X process at the 7 TeV LHC 
in the p T ' and the J/ip rapidity intervals by taking the Set 1 LDME parameters. 
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Table 3: The cross sections for the pp — > J/tp + W — > J/ip + \x + v^ + X process at the 7 TeV LHC 
in the p T ' and the J/ip rapidity intervals by taking the Set 2 LDME parameters. 



set of LDME parameters. In each table, we list the differential cross sections for fi r 



f l f 



in 



J/i> 



and fi r = fif = mw, respectively. We adopt the event select criteria with the cuts of pit > 25 GeV, 
1 77^1 < 2.4, p™ lss > 20 GeV and ra^ > 40 GeV as requested by ATLAS experiment where muon 
and neutrino are the W decay products. In these two tables we list the J/ip production rates in the 
following pW intervals of 8.5 ~ 10 GeV, 10 ~ 14 GeV, 14 ~ 18 GeV, 18 ~ 22 GeV, 22 ~ 30 GeV, 
and the J/ip rapidity intervals of \yj/^\ < 1.0 and 1.0 < \yji$\ < 2.1, separately. The cross section for 
the pp — > J/ip + W — > J/ip + fi + v^ + X process is obtained via the cross section for pp — > J/ip + W + X 
multiplying the branching fraction for the W — > fiu^ decay, which is taken as 10.57% [9]. 

In Figsj3]^a,b) and FigsH)^a,b), we present the the LO and NLO QCD corrected distributions of the 
transverse momentum of J/ip for the pp — > J/ip + W — > J/ip + fj, + z/„ + X process at the 7 TeV LHC 
in the region restricted by the conditions shown in (|0.2p by adopting the Set 1 LDME parameters. 
The renormalization/factorization scales are set to be \i = m T and \i = mw in Figsj3] and FigsJU 
separately. In Figsj3]^a,b) and FigsJ3Ja,b), we give the predictions for the transverse momentum 
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Figure 3: The LO and NLO QCD corrected distributions of the transverse momentum of J/ip for the 
pp —r J/ip + W — * J/ip + /x + v^ + X process at the 7 TeV LHC in the region with the constraints 

fij = rfij, . (aj 



on final particles shown in (|0.2p by taking the Set 1 LDME parameters and [i = fi r 
\yj^\ < 1.0, (b) 1.0 < \yj^\ < 2.1. 
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Figure 4: The LO and NLO QCD corrected distributions of the transverse momentum of J/ip for the 
pp —r J/ip + W—t J/ip + n + Up + X process at the 7 TeV LHC in the region with the constraints 
on final particles shown in (|0.2p by taking the Set 1 LDME parameters and /j, = /j, r = (if = my/, (a) 
\y m \ < 1.0, (b) 1.0 < \yj^\ < 2.1. 
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Figure 5: The LO and NLO QCD corrected distributions of the transverse momentum of J/ip for the 
pp —> J/ip + W — > J/ip + jjl + Un + X process at the 7 TeV LHC in the region with the constraints 



on final particles shown in (|0.2p by taking the Set 2 LDME parameters and [i = \i r = fif 
\y m \ < 1.0, (b) 1.0 < \y J/i} \ < 2.1. 



m 



J/1> 



(a) 



distributions of J/ip in two rapidity intervals on the J/ip: \yj/ip\ < 1-0 and 1.0 < \yj/f\ < 2.1. As 
a comparison, we also depict the contributions from the 1 Sq and 3 S\ and 3 PJ Fock states in 
Figsl3ja,b) and Figs|Ha,b). The short distance part of the differential cross section contributed by 
the 3 PJ 8) (J = 0,1,2) Fock states is negative, and the LDME < O J /^[ 3 P (8) ] > with Set 1 LDME 
parameters is also negative. We can see from Eq. fjO.lH that the differential cross section contributed 
by the 3 Pj (J = 0, 1, 2) Fock states should be positive. These figures show clearly that the differential 
cross section at the LO is significantly enhanced by the QCD corrections. 

In Figsj5]^a,b) and Figsl6^a,b), we give the LO and NLO QCD corrected distributions of the 
transverse momentum of J/ip for the pp — > J/ip + W — > J/ip + \i + v M + X process at the 7 TeV 
LHC taking the Set 2 LDME parameters in the region restricted by the fiducial cuts on final particles 
as shown in f|0.2|) . The renormalization/factorization scales are set to be \i = m T and fi = m\y in 
Figsj5]^a,b) and Figsj6|a,b), separately. We also give the predictions for the transverse momentum 
distributions of J/ip in two J/ip rapidity intervals: \y,j/^\ < 1.0 and 1.0 < \yj/^\ < 2.1. Figj5]^a), 
FigJUJa) are for \yj/^\ < 1.0, and Figj5]^b), Figj6jb) for 1.0 < \yj/^,\ < 2.1. The short distance 
contribution part from the 3 P) ( J = 0, 1, 2) Fock states is negative, while the LDME < J /^[ 3 P O ] > 
with Set 2 LDME parameters is also positive, thus the contribution from 3 Pj (J = 0,1,2) Fock 
states is negative. The negative contributions of 3 Pj (J = 0, 1,2) Fock states are also presented in 
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Figure 6: The LO, NLO QCD corrected distributions of the transverse momentum of J/ip for the 
pp — > J/ifi + W — > J/t/j + /x + v^ + X process at the 7 TeV LHC in the region with the constraints 
on final particles shown in fjQ.2|) by taking the Set 2 LDME parameters and \x = /J, r = [if = rayy. (a) 
\y J/f \ < 1.0, (b) 1.0 < \y JH> \ < 2.1. 

Figsl5ja,b) and Figsj6|a,b) for comparison. 

In summary, we provide the predictions on the J/ip + W associated production within the factor- 
ization formalism of the nonrelativistic QCD at the 7 TeV LHC using the parameters and kinematic 
constraints suggested by the ATLAS experiment. We show the distributions of the transverse momen- 
tum of J/ip for the pp — > J/ip + W — > J/ip + /i + v M + X process in a region constrained by putting 
the cuts on final particles shown in (|0.2p , and using the two set of LDME parameters and two typical 
energy scales separately. Our results show that the theoretical uncertainties come from the choices of 
the scales and the CO LDME parameters. 
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